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We derive transport equations for the propagation of water wave action in the pres- 
ence of a static, spatially random surface drift. Using the Wigner distribution W(x, k, t) 
to represent the envelope of the wave amplitude at position x contained in waves with 
wavevector k, we describe surface wave transport over static flows consisting of two length 
scales; one varying smoothly on the wavelength scale, the other varying on a scale com- 
parable to the wavelength. The spatially rapidly varying but weak surface flows augment 
the characteristic equations with scattering terms that are explicit functions of the cor- 
relations of the random surface currents. These scattering terms depend parametrically 
on the magnitudes and directions of the smoothly varying drift and are shown to give 
rise to a Doppler coupled scattering mechanism. The Doppler interaction in the presence 
of slowly varying drift modifies the scattering processes and provides a mechanism for 
coupling long wavelengths with short wavelengths. Conservation of wave action (CWA), 
typically derived for slowly varying drift, is extended to systems with rapidly varying 
flow. At yet larger propagation distances, we derive from the transport equations, an 
equation for wave energy diffusion. The associated diffusion constant is also expressed 
in terms of the surface flow correlations. Our results provide a formal set of equations 
to analyse transport of surface wave action, intensity, energy, and wave scattering as a 
function of the slowly varying drifts and the correlation functions of the random, highly 
oscillatory surface flows. 



1. Introduction 

Water wave dynamics are altered by interactions with spatially varying surface flows. 
The surface flows modify the free surface boundary conditions that determine the dis- 
persion for propagating water waves. The effect of smoothly varying (compared to the 



wavelength) currents have been analysed using ray theory ( Peregrine (1976) , Jonsson 



fc Stewart (1961), Mei (1979), White (1999), Whitham (1974) and references within). 



(1990)1) and the p rinciple of conservation of wave action (CWA) (cf. Longuct-Higgins 



These studies have largely focussed on long ocean gravity waves propagating over even 
larger scale spatially varying drifts. Water waves can also scatter from regions of un- 
derlying vorticity regions smaller than the wavelength Fabrikant fc Raevsky (1994)| and 
Cerda & Lund (1993). Boundary conditions that vary on capillary length scales, as well 
as wave interactions with structures comparable to or smaller than the wavelength can 
also lead to wave scattering (Chou, Lucas & Stone (1995) , pou, Messiter fc Schultz 
(1993)), attenuation (Chou & Nelson (1994), Lee et al. (1993)), and Bragg reflections 



( phou (1998) , Naciri fc Mei (1988)[ ). Nonetheless, water wave propagation over random 
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static underlying currents that vary on both large and small length scales, and their 
interactions, have received relatively less attention. 

In this paper, we will only consider static irrotational currents, but derive the transport 
equations for surface waves in the presence of underlying flows that vary on both long 
and short (on the order of the wavelength) length scales. Rather than computing wave 
scattering from specific static flow configurations ( perber (1993) , Trulson & Mei (1993) 



Fabrikant & Raevsky (1994)), we take a statistical approach by considering ensemble 



averages over realisations of the static randomness. Different statistical approaches have 
been applied to wave propagation over a random depth (Elter & Molyneux (1972)), third 
sound localization in superfluid Helium films (Kleinert (1990)), and wave diffusion in the 
presence of turbulent flows (Howe (1973), Rayevskiy (1983), Fannjiang & Ryzhik (1999)). 

In the next section we derive the linearised capillary-gravity wave equations to low- 
est order in the irrotational surface flow. The fluid mechanical boundary conditions are 
reduced to two partial differential equations that couple the surface height to velocity 



potential at the free surface. We treat only the "high frequency" limit (Ryzhik, Papani- 



colaou, & Keller (1996)) where wavelengths are much smaller than wave propagation dis- 



tances under consideration. In Section 3, we introduce the Wigner distribution W(x, k, t) 
which represents the wave energy density and allows us to treat surface currents that 
vary simultaneously on two separated length scales. The dynamical equations developed 
in section 2 are then written in terms of an evolution equation for W. Upon expanding 
W in powers of wavelength/propagation distance, we obtain transport equations. 

In Section 4, wc present our main mathematical result, equation (4.1), an equation 
describing the transport of surface wave action. Appendix A gives details of some of 
the derivation. The transport equation includes advection by the slowly varying drift, 
plus scattering terms that are functions of the correlations of the rapidly varying drift, 
representing water wave scattering. Upon simultaneously treating both smoothly varying 
and rapidly varying flows using a two-scale expansion, we find that scattering from rapidly 
varying flows depends parametrically on the smoothly varying flows. In the Results and 
Discussion, we discuss the regimes of validity, consider specific forms for the correlation 
functions, and detail the conditions for doppler coupling. CWA is extended to include 
rapidly varying drift provided that the correlations of the drift satisfy certain constraints. 
We also physically motivate the reason for considering two scales for the underlying 
drift. In the limit of still larger propagation distances, after multipl e wave scatter ing, 
wave propagation leaves the transport regime and becomes diffusive ( Sheng (1995) ). A 
diffusion equation for water wave energy is also given, with an outline of its derivation 
given in Appendix B. 



2. Surface wave equations 

Assume an underlying flow V(x, z) = (J7i(x, z), [/2(x, z), [/^(x, z)) = (U(x, z), [/^(x, z)), 
where the 1,2 components denote the two-dimensional in-plane directions. This static 
flow may be generated by external, time independent sources such as wind or inter- 
nal flows beneath the water surface. The surface deformation due to V(x, z) is denoted 
77(x) where x = (a;,y) is the two-dimensional in-plane position vector. An additional 
variation in height due to the velocity v(x, z) associated with surface waves is denoted 
r]{x,t). When all flows are irrotational, we can define their associated velocity poten- 
tials V(x, z) = (Vx + z9z)<I>(x, z) and v(x, z, t) = (Vx + z9z)<y9(x, z, t). Incomprcssibility 
requires 

A(^(x, z, t) + 5,V(x, z, t) = A$(x, z) + 92$(x, z) = 0, (2.1) 
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Figure 1. The relevant scales in water wave transport. Initially, the system size, observation 
point, and length scale of the slowly varying drift is 0{L), with surface wave wavelength and 
scale of the random surface current of 0(1). Upon rescaling, the system size becomes 0(1), while 
the wavelength and random flow variations are 0{e). 



where A = is the two-dimensional Laplacian. The kinematic condition ( Whitham 
(1974)1 ) applied at z = ^(x) + ?7(x, t) = C(x, t) is 

9*77(x, t) + U(x, C) • VxC(x, t) = C/,(x, z = C) + dM^, ^ - C, t). (2.2) 



Upon expanding (2.2) to linear order in 77 and tp about the static free surface, the right 
hand side becomes 



[/,(x,C) + a,(^(x,C,i) = i7,(x,fy) + 77(x,t)9,[/,(x,77) + a,(^(x,77,t) + 0(772). 



(2.3) 



At the static surface fj, U(x, 77) • Vx^(x) — Uz{x, fj). Now assume that the underlying flow 
is weak enough such that ?7z(x, z w 0) and fj are both small. A rigid surface approximation 
is appropriate for small Froude numbers W^/c^ ~ |Vx?7p ~ J7z(x, 0)/|U(x, 0)| <C 1 (c^ 
is the surface wave phase velocity) when the free surface boundary conditions can be 
approximately evaluated at z = (Fabrikant & Raevsky (1994)). Although we have 
assumed J7z(x, z « 0) = 9z*&(x, z « 0) « and a vanishing static surface deformation 
7?(x) « 0, Vx • U(x, 0) = -9z?7z(x, 0) ^ 0. 

Combining the above approximations with the dynamic bounda ry conditions (derived 
from balance of normal surface stresses at z = ( Whitham (1974) )), we have the pair of 
coupled equations 



i9t77(x, t) + Vx • (U(x, z = 0)7y(x, t)) = lim ^^(^(x, z, t) 



0- 



lim [pdttp{x, z, t) 

z— »0~ 



pU(x, z) • Vx<rf3(x, z, t)] aAi]{x, t) - pgri{y., t) 



(2.4) 



where a and g are the air-water surface tension and gravitational acceleration, respec- 
tively. Although it is straightforward to expand to higher orders in /^(x) and 77(x, i), or to 
include underlying vorticity, we will limit our study to equations (^^) in order to make 
the development of the transport equations more transparent. 



The typical system size, or distance of wave propagation shown in Fig. ^ is of 0{L) with 
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L » 1. Wavelengths however, are of 0{1). To implement our high frequency (Ryzhik 



Papanicolaou, & Keller (1996)) asymptotic analyses, we rescale the system such that all 
distances are measured in units oi L = e~^. We eventually take the limit e — s- as an 
approximation for small, finite e. Surface velocities, potentials, and height displacements 
are now functions of the new variables x — s- x/e, z z/e and t — > t/s. We shall further 
nondimensionalise all distances in terms of the capillary length £c — y/a/gp. Time, 
velocity potentials, and velocities are dimensionalised in units of \/T7[g, \/ gl\., and \fgi'c 
respectively, e.g. for water, [7 = 1 corresponds to a surface drift velocity of ^ 16.3cm/s. 

Since ?72(x, z « 0) w 0, we define the flow at the surface by 

U(x, z = 0) = U(x) + ^/i5U(x/e). (2.5) 

In these rescaled coordinates, U(x) denotes surface flows varying on length scales of 
0(1) much greater than a typical wavelength, while (5U(x/e) varies over lengths of 0[e) 
comparable to a typical wavelength. The amplitude of the slowly varying flow U(x) is 
0(e°), while that of the rapidly varying flow 5U(x/e), is assumed to be of 0{y/e). A more 
detailed discussion of the physical motivation for considering the ^/e scaling is deferred 
to the Results and Discussion. After rescaling, the boundary conditions ( |2.4[ ) evaluated 
at 2: = become 

5t7;(x,i)+Vx- [(U(x) + ^/i<5U(x/£)) 77(x)] = lim a,(^(x,0) 

dt(p{x, t) + U(x) • Vx</?(x, t) + x/i5U(x/e) • Vx(/3(x, t) = eA?7(x, t) - e^\{x, t). 

(2.6) 

Although drift that varies slowly along one wavelength can be treated with characteristics 
and WKB theory, random flows varying on the wavelength scale require a statistical 
approach. Without loss of generality, we choose (5U to have zero mean and an isotropic 
two-point correlation function {SUi{x.)SUj{x')) = i?ij(|x — x'|), where — (1,2) and 
(. . .) denotes an ensemble average over realisations of (5U(x). 
We now define the spatial Fourier decompositions for the dynamical wave variables 

^{^,-h^z^C,t)= J y(q,t)e-'q- '°'^^^^^||'^+^\ 77(x,i)= J r;(q, t)e-''i •^ 



the static surface flows 



(2.7) 

U(x)= J U(q)e-*i '', '5U(^) = J (5U(q)e-''^•''/^ (2.8) 



and the correlations 



i?y (x) = f R,j{q)e-'^^-^, (2.9) 

where q = (<Zi,92) is an in-plane two dimensional wavevector, q = |q| = \/qf + and 
/ = {2tt)~'^ J dqidq2. The Fourier integrals for r/ exclude q = due to the incompress- 
ibility constraint rj{yi, t) = 0, while the q = mode for (p gives an irrelevant constant 



shift to the velocity potential. Note that ip in (2.7) manifestly satisfies (2.1). Substituting 
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(2.8) into the boundary conditions (2.4), we obtain, 

dtrj{\i,t) — i / ry(k — q)U(q) • k — i-^e / 77(k — q/e)(5U(q) • k = i^(k, f) fc tanhe/cft. 

./ q ./ q 

t)^i f U(q) • (k - q) ^(k -Ci)-iV~£ I ^U(q) • (k - q/e) ^(k - q/e) 

./ q ./ q 



'q 

U2 I ^-IN 



= -(efc2+£-l)7y(k). 

(2.10) 

where the (5U(q) are correlated according to 

{5U,{v)5U,{ci)) = i?.,(|p|)<5(p + q). (2.11) 

Since the correlation i?y (x) is symmetric in i j, and depends only upon the magnitude 
|x|, Rij{\p - q|) is real. 



In the case where SJJ = and U(x) = Uq is strictly uniform, equations ( 2.10 ) can 
be simplified by assuming a e~*'^* dependence for all dynamical variables. Uniform drift 
yields the familiar capillary-gravity wave dispersion relation 



u;(k) = ^(fc3 + k) tanh fc/i + Uq • k = n{k) + Uq • k. (2.12) 

However, for what follows, we wish to derive transport equations for surface waves (action, 
energy, intensity) in the presence of a spatially varying drift containing two length scales: 
U = U(x) + V?5U(x/e). 



3. The Wigner distribution and asymptotic analyses 

The intensity of the dynamical wave variables can be represented by the product of two 
Green functions evaluated at positions x±£y/2. The difference in their evaluation points, 
ey, resolves the waves of wavevector |k| ~ 27r/(e?/). Elter & Molyneux (1972) used this 



representation to study shallow water wave propagation over a random bottom. However, 
for the arbitrary depth surface wave problem, where the Green function is not simple, 
and where two length scales are treated, it i s convenient to use th e Fourier representation 
of the Wigner d istribution ( Wigner (1932) , Gerard et al. (1997) , Ryzhik, Papanicolaou 
k Keller (1996)| ). 

Define tp — {ipi, ip2) = (?7(x), z = 0)) and the Wigner distribution: 



W,,(x, k, t) EE (27r)-2 J e^k-y^^ (^^ _ El^t) (x + ^, i) dy 



(3.1) 



where x is a central field point from which we consider two neighbouring points x ± ^ , 
and their intervening wave field. Fourier transforming the x variable using the definition 
(2.7) we find. 



W,,{pM = (27re)-V. ( ? ~ ) V^* 



(3.2) 



The total wave energy, comprising gravitational, kinetic, and surface tension contribu- 
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£^ [|Vx^(x)p + h(x)H +\j^f^ rf^|U(x, z) + z[/,(x, z) + v(x, z)| 



dz\\5{^, z) + zUz{x, z) 



(3.3) 



= - I (fc2 + l)|7y(k)|2 + /jtanhA:/i|(^(k,z = 0)p. 

The energy above has been expanded to an order in 7y(x, t) and </?(x, z, t) consistent with 
the approximations used to derive (2.4). In arriving at the last equahty in (3.3), we have 
integrated by parts, used the Fourier decompositions (2.7) and imposed an impe netrable 
bottom condition at z = —h. The wave energy density carried by wavevector k is (Gerard 



et al. (1997)1 ) 



£:(k,<) = iTr[A(k)W(k,i)] 



(3.4) 



where ^ii(k) = + A22(k) — fc tanh fcft,, v4i2 = ^21 = 0. Thus, the Wigner distribution 
epitomises the local surface wave energy density. 

In the presence of slowly varying drift, we identify W(x,k, t) as the local Wigner 
distribution at position x representing waves of wavevector k. The time evolution of its 
Fourier transform W(p, k, t), can be derived by considering time evolution of the vector 
field if) implied by the boundary conditions (2.4): 



V'i(k,t) +iLj-£(k)V'£(k,t) 



U(q).(k-q5,2)^j(k-q,t) 



(3.5) 



+zV^/q(5U(q) • (k-q5,2/e)Vj(k-q/e,0, 
where the operator L(k) is defined by 

i|k|tanhe|k|/i 



L(k) 



-i(efc2+£-i) 



(3.6) 







We have redefined the physical wavenumber to be k/e so that k ~ 0(1). Upon using 
( ^.5[ ) and the definition (3^), (see Appendix A) 



W',j(p,k,t) =iWa(p,k,i)4, 



- iL, 



£ 2 



+i 



I U(q) . (^-^ + I - q<S,2) W,,{p - q, k + eq/2. 



U(q) 



t) 



q<5j2 Wij (p - q, k - eq/2, t) 



(3.7) 



+iV^JjV{q) . (^-^ + I - W,, (p - q/e, k + q/2, t) 

-iV^ljUiq) . (^-^ - I + ^S,2^ W,,(j> - q/e, k - q/2, t), 

where only the index ^ = 1, 2 has been summed over. If we now assume that W(x, k, t) 
can be expanded in functions that vary independently at the two relevant length scales, 
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functions of the field p (dual to x) can be replaced by functions of a slow variation in p 
and a fast oscillation ^/e; p ^ p + ^/e- 

This amounts to the Fourier equivalent of a two-scale expansion in which x is replaced 
by X and y = x/e (Elyzhik, Papanicolaou, & Keller (1996)). The two new independent 
wavevectors p and ^ are both of 0(1). Expanding the Wigner distribution in powers of 
^/e and using p ^ p + ^/e, 

W(p, k, t) -> Wo(p, I, k, t) + ^/iWl/2(p, k, t) + £Wi(p, I, k, t) + 0(£3/2), (3.8) 

we expand each quantity appearing in (^^) in powers of ^/e and equate like powers. 
Upon expanding the off-diagonal operator L(— k/e + p/2) = £^^Lo(k) +Li(k, p) 
where 



Lo(k) = 








and 



f{k) 



Li(k,p) 



hk + sinh kh cosh kh 
2k cosh^ kh 




(3.10) 



3.1. Order e 



terms 



The terms of 0{e ^) in (3.7) are 

Wo(p,e,k,i)LS(k+)-Lo(k_)Wo(p,^,k,i) =0, k± = k±| (3.11) 

To solve (^Tll), we use the eigenvalues and normalised eigenvectors for Lq and its complex 
adjoint lJ, 

/ zrv/^(k)72 \ 
Tf2(k) — ij,hr = 2 ' TQCk) + i-f, - 

\ VMk) / 

where a(k) = /^^^"^ , r = ±1, and 17 —> is a small imaginary term. A Wo(p, |, k, t) 
fc^ -I- 1 

that manifestly satisfies (3.12) can be constructed by expanding in the basis of 2 x 2 
matrices composed from the eigenvectors: 



VMk) 

y^(k)72 



(3.12) 



Wo(p,^,k,i) = «r.'(P,k,<)b,(k_)bt,(k+) 



r,r' — ± 



(3.13) 



Right [left] multiplying ( 3.11 ) (using ( p.l3| )) by the eigenvectors of the adjoint prob- 
lem, CT-(k_) [c^(k+)], we find that = a |- = 0, and a (x, k, i) = a_(x, k, i) = 

a++(x, — k, f) = a-|-(x, — k, f). Furthermore, a+,a_ ^ only if ^ = 0. Thus Wq has the 
form 

Wo(p, I, k, t) = Wo(p, k, t)m- (3.14) 
From the definition of Wq, we see that the (1,1) component of Wq is the local envelop 
of the ensemble averaged wave intensity |?7(x, k, t)p ~ a+(x, k, t)a(k). Similarly, from the 
energy (Eq. (p]j)), we see immediately that the local ensemble averaged energy density 



(f(x,k,t)) = Aii(k)a(k)(a(x,k,i)) +.l22(k)(a(x,k,i)) 
= fi(k)(a(x,k,i)). 



(3.15) 
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Therefore, since the starting dynamical equations are hnear, we can identify (a(x, k, t)) as 
the ensemble averaged local wave action associated with waves of wavevector k ( Henyey 
et al. (1988)). The wave action (a(x, k, t)), rather than the energy density (£'(x, k, i)) 
is the conserved quantity (Longuct-Higgins & Stewart (1961), Mei (1979), Whitham 
(1974)D. 



The physical origin of 7 arises from causality, but can also be explicitly derived from 
considerations of an infinitesimally small viscous dissipation (Chou, Lucas & Stone 
(1995)). Although we have assumed 7 — > 0, for our model to be valid, the viscosity 
need only be small enough such that surface waves are not attenuated before they have 
a chance to multiply scatter and enter the transport or diffusion regimes under consid- 
eration. This constraint can be quantified by noting that in the frequency domain, wave 



dissipation is given by 7 = 2vk (Landau (1985)) where v is the kinematic viscosity and 



Cg{k) = |Vkr!(fc)| 



(3.16) 



is the group velocity. The corresponding decay length fc^ ^ ~ Cg(fc) / [vk"^) must be greater 
than the relevant wave propagation distance. Therefore, we require 



^^M^»(l,e-) 



(3.17) 



for (transport, diffusion) theories to be valid. The inequality ( 3.17 ) gives an upper bound 
for the viscosity 



vk^ < {ecg{k/e),e^cg{k/e)) 



(3.18) 



which is most easily satisfied in the shallow water wave regime for transport. Otherwise 
we must at least require v < o{y^). The upper bounds for v (and hence 7) given above 
provide one criterion for the validity of transport theory. 



3.2. Order e terms 
Collecting terms in (3.7) of order e~^^^, we obtain 

Wi/2(p, ^, k, i)LS(k+) - Lo(k_)Wi/2(p, I, k,t)+ f U(q) • |Wi/2(p - q, I - q, k, t) 

- f (5U(q) • k_ Wo(p, ^ - q, k + q/2, t) + f JU(q) • k+ Wo(p, | - q, k - q/2, t) 

- / (5U(q) • q [Wo(p, I - q, k + q/2, t)S + SWo(p, ^ - q, k - q/2, t)] = 

Jq 

(3.19) 

where S = q 1 

Similarly decomposing W1/2 in the basis matrices composed of bT-(k_)b|,(k_|_) (as in 



3.13), substituting Wo(p, 0, k, t)(5(^) from (3.13 into (3.19), and inverse Fourier trans- 
forming in the slow variable p, we obtain 

w , k r t) V ^^^^^ • r.,.'(x,l,k,t)b.(k_)bUk+) 
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where 



Tr,r' (x, ^ , k, t) EE k_a^- (x, k+ , t)ct (k_ )b^- (k+) - k+ar (x, k_ , t)hl (k_ )c^/ (k+) 

+ 1 J2 k+, i)ct (k_)b^(k+) + a^x, k_, t)bt (k_)c,,(k+)] 



3.3. Order terms 
The terms of order e° in (pj) read 



(3.21) 



Wo(p, k, t) = zWo(p, k, i)Ll (-p) - *Li(p)Wo(p, k,t)~z k- U(q)q • VkWo(p - q, t K t) 

Jq 

+t I U(q).pWo(p-q,ek,t)-* / U(q).q[SWo(p-q,|,k,t)+Wo(p-q,C,k,t)S] 
+i [ <5U(q)-k+Wi/2(p,^-q,k-q/2,t)-z / (5U(q) • k_ Wi/2(p, ^ - q, k + q/2, t) 

^q J c{ 

^ I <5U(q).q[SWi/2(p,^-q,k + q/2,t)+Wi/2(p,|-q,k-q/2,t)S] 

Jq 

+^WlLt - zLoWl + f U(q).|Wi(p-q,^,k,t). 

Jq 

(3.22) 

To obtain an equation for the statistical ensemble average (a+(x, k, i)), we multiply 
(B.22) by c^(k) on the left and by c+(k) on the right and substitute W1/2 from equation 
(B.20). We obtain a closed equation for a(x, k, i) = (a+(x, k, t)) (we henceforth suppress 
the (. . .) notation for a(x, k, t) and £(x, k, t)) by truncating terms containing Wi. Clearly, 
from ( [3.12| ), c^(k)(iWiLj — iLoWi)c_|_(k) — 0. Furthermore, we assume (^Wi(p — 
q, ^, k, t)) « which follows from ergodicity of dynamical systems, and has been used in 



the propagation of waves in random media (see Ryzhik, Papanicolaou, & Keller (1996) 



Bal ct al. (1999)). The transport equations resulting from this truncation are rigorously 



justified in the scalar case (3pohn (1977), Erdos & Yau (1998)). 



4. The surface wave transport equation 

The main mathematical result of this paper, an evolution equation for the ensemble 
averaged wave action a(x, k, t) follows from equation (3.22) above (cf. Appendix A) and 
reads, 

d(x, k, t) + Vkij-'(x, k) • Vxa(x, k, t) — Vxw(x, k) • Vka(x, k, t) 



-I](k)a(x, k, t) + / a(q, k)a(x, q, t), 

J a 



(4.1) 



where 



w(x, k) = yj{k^ + k) tanhfc/i + U(x) • k = Vl{k) + U(x) • k. (4.2) 

The left hand side in (4.1) corresponds to wave action propagation in the absence of 
random fluctuations. It is equivalent to the equations obtained by the ray theory, or 
a WKB expansion (see section 5.1). The two terms on the right hand side of (4.1) 
represent refraction, or "scattering" of wave action out of and into waves with wavevector 
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k respectively. In deriving ( [4.lD we have inverse Fourier transformed back to the slow field 
point variable x, and used the relation (Q!(k) — /(/c)q!~^ (k))k = Vki^(k). To obtain ( |4.lD , 
we assumed Rij{q)qi = Rij{q)qj = 0, which would always be valid for divergence-free 
flows in two dimensions. Although the perturbation S\J is not divergence-free in general, 
V-JU(x, z = 0) = —dz5Uz{x.,0) 7^ 0, using symmetry considerations, we will show in 
section 5.2 that Rij{q)qi — Rij{(\)qj — 0. 
Explicitly, the scattering rates are 

E(k) = 2^ / g,i?y(q-k)fc, b^k)c,(q)bt(q)c+(k)<5(™(x,Tq) -^(x,k)) 
cr(q,k) = 27r ^ T(?ji?ij(Tq - k)A:j|b]:(Tq)c+(k)p(5 (™(x,q) - c^(x,k)) 



where 



, t lr.\\.Ur.\r. (Ta(k) + Q;(q))(Ta(q) + a(k)) 
bUk)c.(q)bt(q)c+(k) = A^^iy^) 

Ibt(k)c.,(q)p + 



(4.3) 



(4.4) 



4a(k)a(q) 

Physically, S(k) is a decay rate arising from scattering of action out of wavevector k. 
The kernel CT(q, k) represents scattering of action from wavevector q into action with 
wavevector k. Note that the slowly varying drift U(x) enters parametrically in the scat- 
tering via a;(x, k) in the ^— function supports. The arguments ti^(x, k) in the 5— functions 
mean that we can consider the transport of waves of each fixed frequency ujq = ti;(x, k) 
independently. 

The typical distance travelled by a wave before it is significantly redirected is defined 
by the mean free path 

The mean free path described here carries a different interpretation from that considered 



in weakly nonlinear, or multiple scattering theories ( Zakharov, L'vov fc Falkovich (1992) ) 
where one treats a low density of scatterers. Rather than strong, rare scatterings over 
every distance ^^jp ^ 0(1), we have considered constant, but weak interaction with 
an extended, random flow field. Although here, each scattering is 0(e) and weak, over 
a distance of 0(1), approximately interactions arise, ultimately producing ^j^jp ^ 
0(1). 



5. Results and Discussion 

We have derived transport equations for water wave propagation interacting with 
static, random surface flows containing two explicit length scales. We have further as- 
sumed that the amplitude of 515 scales as with (3 = 1/2: The random flows are 
correspondingly weakened as the high frequency limit is taken. Since scattering strength 
is proportional to the power spectrum of the random flows and is quadratic in 5U, the 
mean free path can be estimated heuristically by ^j^jj^ ~ Cg{k) /Yi{k)£^^'^^ . For fi > 1/2, 
the scattering is too weak and the mean free path diverges. In this limit, waves are nearly 
freely propagating and can be described by the slowly varying flows alone, or WKB the- 
ory. If /3 < 1/2, i^j^p and the scattering becomes so frequent that over a propagation 
distance of 0(1), the large number of scatterings lead to diffusive (cf. Section 5.4) be- 
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haviour (3hcng (1995)). Therefore, only random flows that have the scaling (3 = 1/2 
contribute to the wave transport regime. 

We also note that f3 > precludes any wave localisation phenomena. In a two- 
dimensional random environment, the localisation length over which wave diffusion is 



inhibited is approximately (Sheng (1995)) 



'loc mfp 



exp 



-1-2/3 



exp [e 



-2f3\ 



(5.1) 



As long as the random potential is scaled weaker {(3 > 0), li^^ oo, and strong lo- 
calisation will not take hold. In the following subsections, we systematically discuss the 
salient features of water wave transport contained in Eq. (4.1) and derive wave diffusion 
for propagation distances > 0(1). 

5.1. Slowly varying drift: U(x) ^ 0, S\J = 
First consider the case where surface flows vary only on scales much l arge r than the 



longest wavelength 27r/fc considered, i.e., SU = 0. The left-hand side in (4T) represents 
wave action transport over slowly varying drift and may describe short wavelength modes 
propagating over flows generated by underlying long ocean waves. 



We first demonstrate that the nonscattering terms of the transport equation (4.1) is 
equivalent to the res ults obtained by ray theory (WKB e xpansion) and conservation of 
wave action (CWA) (|Longuet-Higgins fc Stewart (1961)| , |Mei (1979)|, [Peregrine (1976 ' 



White (1999)1 IWhitham (1974)1) . Assume the WKB expansion ( [Keller (1958)| , [Bender fc 
Orszag (1978)[ ) 

?fe = A„ (x, i)e^^^(-'*)/^ and - (x, t)e'^^^-'^/' , 



(5.2) 



with smoothly varying A,j and A^p. Upon using the above ansatz in (3.1) and setting 



0, we have a(x,k,i) = | A|2(x, i)(5(k - Vx5'(x, t)) where \A\'^ = 2a{k)\A^ 



2a ^{k)\Ari\'^. Substitution of this expression for a(x, k, i) into (4.1), we obtain the fol- 
lowing possible equations for S'(x, t) and |Ap(x, t) 



5t|^|2(x,t)+V, 



at5 + c^(x, Vx^) = o, 

(|A|2VkC^(x,Vx^)) -0. 



(5.3) 
(5.4) 



The first equation is the eikonal equation, while the second equation is the wave action 
amplitude equation. Recalling that = Q;(fc)|Ap/2, we obtain the following transport 
equation for the height amplitude: 

|2 



I A, 



^a(Vx5) 



«(Vx5) 



Vkc^(x, Vx5) 



0. 



(5.5) 



Equation (5_J) is the same as Eq. (8) of White (1999) , except that his Vt is replaced here 



with a due to our inclusion of surface tension. 

Wave action conservation can be understood by noting that 

^a(X(i),i^(i),i)=0, 

where the characteristics {X{t),K{t)) satisfy the Hamilton equations 
dX{t) ^ , , dK{t) 



dt 



\7MXit),K{t)), 



and 



dt 



VMX{t),K{t)). 



(5.6) 



(5.7) 



The solutions to the ordinar y differential equations ([5.7] ) are the characteristic curves 
used to solve (fj) and ^) ( [Courant fc Hilbert {196if ). 
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5.2. Correlation functions and conservation laws 



We now consider the case where 5\5 ^ 0. The scattering rates defined by (4.3) depend 
upon the precise form of the random flow correlation i?y . There are actuaUy six additional 
terms in (4.3) in the calculation of <j and S, which vanish because 

2 

^%(q)9j=0 for z = 1,2. 



(5.8) 



We prove relation (5.8) provided that (5[/z(k, fc^) and ^[/^(k, ~kz) have the same proba- 
bility distribution. Thus, 



(^(7,(p,p,)(5C/,(k,fc,)) = ((5[/,(p,p,)^[/,(k,-fc,)) 



(5.9) 



This symmetry condition is reasonable, and is compatible with the divergence-free con- 
dition for i5U in three dimensions. We show that Hypothesis ( ^.9[ ) implies (5.8) by first 
using incompressibility Y^j=i '^t^i(k, kz)kj + (5C/z(k, kz)kz — 0: 



^5(p + k)i?,,(k)fc, = ^((5f/,(p,0)(5C/,(k,0)A:,) 

^((5[/,(p,0) / 5U,{\^,kz)k,) 

..■ — 1 ^ 

{5U,{p,Q)5Uz{\^, kz))kzdkz 



0, 



where the last equality follows from (5_^). Thus, (5^) is verified, and ( |4.3D derived. 
The form Rij(\c^)qi — i?y(|q|)gj — 0, requires the correlation function to be transverse: 



RM)^R{q) 



HI] 



(5.10) 



where R{q) is a scalar function of q. The correlation kernels in the scattering integrals 
can now be written as 



qiRij{\TCi - k|)fc, 



i?(|rq-k|) 



qk 



q • (rq — k)k • (rq — k) 



|Tq-k|2 



^ i^(lrq-kl) ^,^,^.^,^ 



(5.11) 



where 9 denotes the angle between q and k. The scattering must also satisfy the support 
of the (5-functions; for U(x) = only |q| = |k| satisfy the the (5— function constraints. In 
the presence of slowly varying drift, the evolution of a(x, |k| ^ |q|) can "doppler" couple 
to a(x, q, t). 



It is straightforward to show from the explicit expressions (O) that 

I](k)= / a(k,q). 
J a 



(5.12) 



This relation indicates that the scattering operator on the right hand side of (4.1) is 
conservative: Integrating (4.1) over the whole phase space yields 



d_ 
di 



a(x, k, t) — 0. 



(5.13) 
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Figure 2. (a). Contour plot of (^(q). Each grayscale corresponds to a different constant value 
of i-o{q) = u!(k) = iiJo- (b). The band of q that satisfies 0.625 < uo < 0.6625. Wavevectors q and 
k that lie in this band can couple a(x, k, t) to a(x, q, t) via wave scattering. 



Equation ( 5.13 ) is the generalization of CWA to include scattering of action from rapidly 
varying random flows 6XJ{x/e). Although a(x, k, t) is conserved, the total water wave 
energy £{x,'k,t) = r2(k)a(x, k, f) will not be conserved. For example, if [/(x) is small 
enough such that the J— function in the (T(q, k) integral is triggered only when r — +1, 



dt 



d 
dt 



j^jM^. k) - k • U(x)]a(x, k, i) = ^ ^ k • U(x) a(x, k, i) ^ 



(5.14) 



This nonconservation results from the energy that must be supplied in order to sustain 
the stationary underlying flow. For small U (x), the quantity ti'(x, k)a(x, k, t) is conserved. 
In that case, the evolution of cj(x, k)a(x, k, t) obeys an equation identical to (4.1). When 
there is doppler coupling with t = —1, an additional term arises and a;(x, k)a(x, k, t) is 
no longer conserved under scattering. 

5.3. Doppler coupled scattering 

In addition to the correlation functions, the wave action scattering terms involving S(k) 
and integrals over cr(q, k) depend also on the support of the (5— function. Consider action 
contained in water waves of fixed wavevector k. When U(x) = 0, only r = +1 terms 
contribute to the the integration over q as long as |q| = |k|. In this case, we can define 
the angle q • k = fc^ cos 9 and reduce the cross-sections to single angular integrals over 

/ 9 \ sin"^ 9 

g,i?y(|q-k|)q,=i? |2fcsin-| — r = +1. (5.15) 
V ^ / 4 sm I 

In this case (U(x) = 0), assuming i?(|q|) is monotonically decreasing, the most important 
contribution to the scattering occurs when q and k are coUinear. 

When U(x) ^ 0, and r = +1, the sets of q which satisfy ri(q) + U(x) • q = 0(k) + 
U(x) • k = ojo trace out closed ellipse-like curves and are shown in the contour plots 
of i^(q) in Figure ^(a). The parameters used are U(x) • ki = — 0.5fci and h = oo (the 
— ki, — qi directions are defined by the direction of U(x)). 



Each grayscale corresponds to a curve defined by fixed w(k) = luq. All wavevectors q 
in each contour contribute to the integration in the expressions for E(k) and w(q, k). 
Thus, slowly varying drift can induce an indirect doppler coupling between waves with 
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difFerent wavenumbers, with the most drastic coupling occurring at the two far ends of 
a particular oval curve. For example, in Figure ^(b), the dark band denotes q such that 
w(q) = ijjQ when 0.625 < ojq < 0.6625. The wavevectors q w (—0.3, 0) and q sa (0.8, 0) 
are two of many that contribute to the scattering terms. Therefore, the evolution of 
a(x, k K, (— 0.3,0), i) also depends on a(x, q « (0.8,0),i) via the second term on the 
right side of (pi]). 

Provided U(x) is sufficiently large, the r = — 1 terms can also contribute to scatter- 
ing. The dissipative scattering rate S(k)a(x, k, t) will change quantitatively since addi- 
tional q's will contribute to S(k). However, this decay process depends only on k and 
is not coupled to a(x, |q| ^ |k|,i). Wavevectors q that satisfy the 5— hmction in the 
(T(q, k)a(x, q, <) term will, as when r = -1-1, lead to indirect doppler coupling. This oc- 
curs when w(q) = — cjq and, as we shall see, allows doppler coupling of waves with more 
widely varying wavelengths than compared to the r = -1-1 case. Observe that if r = — 1 
terms arise, the drift frame energy a(x, k, t)aj(x, k) is no longer conserved. Figure ^(a) 
plots uj{qi,q2 = 0) for [/(x) = 1 < \/2, C/(x) = %/2, and C/(x) = 1.6 > V2. Since luq 
and a;(q) are identical functions, —loq can take on values below the upper dotted line 
(wo ^ 0.22 for U — 1.6). Therefore, coupling for r — —1 and (72 = occurs for values of 
—Wo between the dotted lines. Note that depending upon the value of wo, coupling can 
occur at two or four different points q ~ ((71, 0). Figure ||(b) shows a contour plot of |w(q)| 
as a function of (qi, (72)- A level set lying between the dotted lines in (a) will slice out two 
bands; one band corresponds to all values of k that couple to q lying in the associated 
second band. The two bands determined by the interval 0.414 < —uoq < 0.468 are shown 
in Fig. ||(c). For any k lying in the inner band of Fig. ^(c), all q lying in the outer band 
will contribute to doppler coupling for t — —1, and vice versa. As —uq is increased, the 
inner(outer) band decreases(increases) in size, with the central band vanishing when —loq 
approaches the upper dotted line in (a) where the r = — 1 coupling evaporates. If —loq 
is decreased, the two bands merge, then disappear as —ujq reaches the lower limit. Fig. 
||(d) is an expanded view of the two bands for small 0.0756 < —luq < 0.1368. Note that 
a small island of q or k appears for very small wavevectors. The water wave scattering 
represented by cr(q, k) can therefore couple very long wavelength modes with very short 
wavelength modes (the two larger bands to the right in Fig. ||(d)). However, the strength 
of this coupling is still determined by the magnitude of qiRij{\ci — k|)fcj, which may be 
small for |q — k| large. 



The depth dependence of doppler coupling will be relevant when hq, hk < I where q 
and k are the magnitudes of the wavevectors of two doppler-coupled waves. For t = +1, 
finite depth reduces the ellipticity of the coupling bands, resulting in weaker doppler 
effects. Since the water wave phase velocity decreases with h, a finite depth will also 
reduce the critical J7(x) required for r = — 1 doppler coupling. For small f/(x), it is clear 
that the (5— functions associated with the t = — 1 terms in cr(q, k) are first triggered when 
the q and k are antiparallel, U • k = — fc|U|, U • q = -|-g|U|. 

Figure |^(a) shows the phase velocity for various depths h. In order for t — —1 to 
contribute to scattering, U ^ c^{k] h). For U ~ 1.6, this condition holds in the h = 00 case 
for 0.5 ^ k <2 (the dashed region of c^(fc, 00)). Recall that our starting equations (2.4) 
are valid only in the small Froude number limit. However, for water waves propagating 
over infinite depth, r = — 1 coupling requires U > Umin = n:iini;{c0(A:)}, with c^(krain) — 
22cm/s. Therefore, in such "supersonic" cases, where r = — 1 is relevant, our treatment 
is accurate only at wavevectors k* such that U ^ c^(k*\ h), e.g., the thick solid portion 
of c</,(fc;oo). For U > Umin, the r = — 1 term can couple wavevectors g « <C kmin 
with k ~ 2 — 3 ^ kmin- The rich t — —1 doppler coupling displayed in Figures || is 
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Figure 3. Conditions for doppler coupling when r — —1. (a). Plot of uj{qi, 52 = 0; ft = 00) for 
U = 1, U = V2, and U = 1.6. Only ior U > V2 does uj{qi,q2 = 0; ft = 00) < 0. (b). Contour 
plot of |aj(q)|. Each grayscale corresponds to a different constant value of CL)(q) = '^(k) = — c^^o. 
(c). The bands of q satisfying 0.414 < —ujo < 0.468. (d). An expanded view of the coupling 
bands for 0.0756 < —ujo < 0.1368. Note that wavenumbers of very small modulus can couple 
with wavenumbers of significantly larger modulus. 



particular to water waves with a dispersion relation a;(q) that behaves as g'^/^jU • q, or 
q^^"^ depending on the wavelength. Doppler coupling in water wave propagation is very 
different from that arising in acoustic wave propagation in an incompressible, randomly 
flowing fluid ( |Howe (1973)| , [Fannjiang fc Ryzhik (1999)| , [Vedantham fc Hunter (1997)D 
where a;(q) = Cs|q|. An additional doppler coupling analogous to the t ~ — 1 coupling 
for water waves arises only for supersonic random flows when C/(x) ^ c^, independent of 
q. In such instances, compressibility effects must also be considered. 

Figure |(b) plots the minimum drift velocity Umin{h) where r = — 1 doppler coupling 
first occurs at any wavevector. The wavevector at which coupling first occurs is also 
shown by the dashed curve. For shallow water, h <C \/3, Umin{h) oc \fh and very long 
wavelengths couple first (small k{Umin))- For depths h > \/3 (~ 3cm for water), the 
minimum drift required quickly increases to U*{oo) = V2, while the initial coupling 
occurs at increasing wavevectors until at infinite depth, where the first wavevector to 
doppler couple approaches k 1 (in water, this corresponds to wavelengths of '--^ 6.3cm). 
The conditions for r = — 1 doppler coupling outlined in Figures ^ and ^ apply to both 
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phase velocity cXk;h) U^^.^ for T=-l coupling 




012345 1 2 3 4 

k h 

Figure A. U > c^(k) is required for t = —1 coupling, (a). The phase velocity C0(fc) for various 
depths h. The velocity shown by the solid horizontal line U « 1.6 > c,j,{k; h = oo) for 0.5 ^ k < 2. 
(b). The minimum Umin(h) required for existence of r = — 1 coupling at any wavevector k, and 
the wavevector k{Umin) at which this first happens. 

E(k) and CT(q, k), with the proviso that q and k are paraUel for S(k) and antiparallel 
for cr(q, k). However, even when U < Umin such that only r — +1 apphes, the set of q 
corresponding to a constant value of a;(k) = tJoi traces out a noncircular curve. There is 
doppler couphng between wavenumbers g 7^ fc as long as C/ ^ 0. 



5.4. Surface wave diffusion 



We now consider the radiative transfer equation (4.1) over propagation distances long 
compared to the mean free path ^jj^^jp- Imposing an additional rescaling and measuring all 

distances in terms of the mean free path, we introduce another scaling e~^, proportional 
to the number of mean free paths travelled. Since f3 = 1/2, transport of wave action 
prevails when 0{e) < |x| ~ 0(1), while diffusion holds when 0(e^^) ~ |x| < i^^^. 



Since waves of each frequency satisfy (4.1) independently, we consider the diffusion of 



waves of constant frequency u>o- To derive the diffusion equation, we assume for simplicity 
that U is constant and small such that ujo + ijj(x, q) ^ (the r = — 1 terms are never 



triggered by the 5— functions). Expanding all quantities in the transport equation (4.1) 
in powers of e, we find 

flo + U- VxCo- Vx-D- Vxflo = 0. (5.16) 
Th e der ivation of this equation is given in Appendix B. The diffusion tensor D is given 



(|B7D: 



i n (|B 12 ) and is a function of the power spectrum Rij. The effective drift U is given by 



- /i,Vkt^(k)J(k-u + r!(k)-^o) 

Up to a change of basis, we can assume that U = Uei, where t7 > 0. Then the set of 
points k • U + r2(k) — = is symmetric with respect to the xi— axis and U is parallel 



to U. Also notice that the total energy given in (3.15) is asymptotically conserved in the 



Water wave transport 



17 



diffusive regime. Indeed, the total energy variations are given by (5.14). Assuming that 
all water waves have frequency ujq, we have in the diffusive regime 



-£ = -- 
dt dt 



i 



J J k-U(x)a(x, k,t) 
y"ao(x,t))^k-U(5(c^o-c^(k)), 



since U is constant. Recasting the diffusion equation as cto — — Vx • (Uao + D • VxOo), 
we deduce that 

J ao(x,t) = 0, 

which conserves the total energy £. 

Now consider the simplified case U = 0, ft, = oo and rioo(k) = \/k^ + k. Since U = 



U = 0, (5.15) holds and we have for all k, 



/ gji?,j(|q- k|)(7jq = 0. 



(5.18) 



We deduce that the corrector x in 



is given by 



^^^^ ^ Vk»oo(k) ^ |Vkr?oo(k)| ^^ 



S(fc) S(fc) " S(fc) 

where k = fck. The isotropic diffusion tensor D is thus given by 



D 



1 



Vq^2oo(q) I ^qq"^ (5(^^00 (q) - ^o) = 



I, 



(5.19) 



m)V^oJq ^ V 2S(fc) 

where I is the 2x2 identity matrix. Thus, the diffusion equation for ao(x, t) assumes the 
standard form ( ^heng (1995) ) 



ao 



cg(fc) 
2I](/c) 



Aoo = 0. 



(5.20) 



6. Summary and Conclusions 

In this paper, we have used the Wigner distribution to derive the transport equations 
for water wave propagation over a spatially random drift composed of a slowly varying 
part U(x), and a rapidly varying part ^/eS\J{x./e). The slowly varying part determines 
the characteristics on which the waves propagate. We recover the standard result obtained 
from WKB theory: conservation of wave action. Provided Rij{q)qj — 0, we extend CWA 
to include wave scattering from correlations Rij of the rapidly varying random flow. 
Evolution equations for the nonconserved wave intensity and energy density can be read- 
ily obtained from (4.1). Moreover, conservation of drift frame energy a(x, k, i)a;(x, k) 
requires small U < Umin and absence of r = — 1 contributions to scattering. 

Explicit expressions for the scattering rates S(k) and cr(q, k) are given in Eqs. ( [4.3D . 
For fixed w(k), we find the set of q such that the (5— functions in (4.3) are supported. 
This set of q indicates the wavevectors of the background surface flow that can me- 
diate doppler coupling of the water waves. Although widely varying wavenumbers can 
doppler couple, supported by the (5— function constraints, particularly for r — —1, the 
correlation i?ij(|q— k|) also decreases for large |q— k|. For long times, multiple weak scat- 
tering nonetheless exchanges action among disparate wavenumbers within the transport 
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regime. Our collective results, including water wave action diffusion, provide a model for 
describing linear ocean wave propagation over random flows of different length scales. 



The scattering terms in (4.1) also provide a means to correlate sea surface wave spectra 



to statistics Rij of finer scale random flows. 



Although many situations arise where the underlying flow is rotational (White (1999)), 
the irrotational approximation used simplifies the treatment and allows a relatively simple 
derivation of the transport and diffusive regimes of water wave propagation. The recent 



extension by White (1999) of CWA to include rotational flows also suggests that an 
explicit consideration of velocity and pressure can be used to generalise the present 
study to include rotational random flows. Other feasible extensions include the analysis 
of a time varying random flow, as well as separating the underlying flows into static and 
wave dynamic components. 

The authors thank A. Balk, M. Moscoso, G. Papanicolaou, L. Ryzhik, and I. Smol- 
yarenko for helpful comments and discussion. GB was supported by AFOSR grant 49620- 
98-1-0211 and NSF grant DMS-9709320. TC was supported by NSF grant DMS-9804780. 



Appendix A. Derivation of the transport equation 

Some of the steps in the derivation of (4.1) are outlined here. By taking the time 
derivative of Wij in ( |3.2D and using the definition (3.5) for ^, we obtain 



(27re)2t^„(p,k,i) = (2^£)2iM/,,(p,k)L| 



k , P 



k p 



^ f U*(q) 
-tV^ I <5U*(q) 

Jq 



e 2 
k p 



{2TT£fiLu 
k P 



k^P 

e 2 



U(q) . + ^ - q5,2 hA. -- + ^ - q -T - ^ 



e 



e 2 

k^P 

e 2 



k P 



k P 



k , p q> 



e 2 e 

k p q 

e 2 e 



k P q 



- + ^ - -fe — + ^ - - w^n — ^ 



Sj2 i^i 



e 2 e 

k^P 

e 2 



k P 



e 2 
(Al) 

To rewrite the above expression as a function of Wij only, we relabel appropriately, 



e.g., 



e 2 



k^ p^ 



k ^ P k' p' 



(A 2) 



for the third term on the right hand side of (Al). Similarly relabelling for all relevant 
terms yields the integral equation (pT7|). 

The 0(e~"'^/^) terms of ( ^.7[ ) determine Wi/2- Decomposing 

(A3) 



Wi/2(p,^,k)= 4!r''(P,l,k)b,(k_)b^,(k+) 



r,r' — ± 



and substituting into ( ^.igj ) we find the coefficients a^^^?^ , where in this case a^^^^ , aL^{.^^ 7^ 
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0. Due to the nonlocal nature of the third term on the right of (3.19), we must first inverse 
Fourier transform the slow wavevector variable back to x. 

To extract the 0(e°) terms from ( |3.7D we need to expand L to order e°, the Li term. 
Similarly, the terms W(p — q, ^, k ± eq/2) must be expanded: 

W(p - q, I , k ± eq/2) = W(p - q, I, k) ± |q • VkW(p - q, k) + 0{e'). (A 4) 

The eq • VkW(p — q, ^, k) terms combine with the — £~^U(q) • k_ + e^^U(q) • k-|_ terms 
from the third and fourth terms in ( |3.7| ) to give the third term on the right of Eq. (3.22). 
The (5U-dependent, order e° terms (the sixth, seventh, and eighth terms on the right side 
of ( 3.22D ) come from collecting 



± VISU{q) . ( ± 1 ) ViWi/2(p, I - q, k ± q/2) 



(A 5) 



from the last two terms in (3.7). The ensemble averaged time evolution of the Wigner 
amplitude ao-(x, k) can be succinctly written in the form: 



a+(x, k, t) — Vxti'(x, k) • Vka+(x, k, t) + Vkw(x, k) • Vxa+(x, k, t) 



S+,^(k)a^(x,k,t) + / (T+^,,(q,k)a^(x,k, 



t). 



(A 6) 

Using the form for Wq found from (3.11) in ( |3.19 ) to find W1/2, we substitute into 
( ^.22 ) to find (4.1), the transport equation for one of the diagonal intensities of the Wigner 
distribution. We have explicitly used eigenbasis orthonormality b]^(k) •Cr'(k) — 5t,t' and 
the fact that a_(x, k, <) = a4-(x, — k, t). 



Appendix B. Derivation of the diffusion equation 

The derivation of diffusion of water wave action is outlined below and follows the 



establi shed mathematical treatment of Larsen & Keller (1974) and Dautray & Lions 



1993)| . For simplicity we assume that the flow U is constant and small enough so that 
for a considered range of frequencies, the relation uj{q) + t^(k) = is never satisfied 
for any k and q 9^ 0. The diffusion approximation is valid after long times and large 
distances of propagation X(see Fig. |l|) such that the wave has multiply scattered and its 
dynamics are determined by a random walk. We therefore rescale time and space as 

~ t ^ X X 

t^^, x=-. Bl 

The small parameter e in this further rescaling represents the transport mean free path 
and not the wavelength as in the initial rescaling used to derive the transport 
equation. We drop the tilde symbol for convenience and rewrite the transport equation 
in the new variables: 

1 If 
a,(x,k,i) + -Vkt^(k)-Vxae(x.k,t) = - / Q(q, k)(ae(x, q, i)-ae(x, k, i))(5(tj(q)-tj(k)), 

(B2) 

with obvious notation for Q(q, k). Since the frequency is fixed, the equation is posed 



for k satisfying w(k) = ujq. The transport equation assumes the form (B2) because the 
scattering operator is conservative. Since U ^ 0, wave action is transported by the flow, 
and diffusion takes place on top of advection. Therefore, we introduce the main drift U, 
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which will be computed explicitly later, and define the drift-free unknown ac(x, k, t) as 

a,(x,k,i) = a,(x + — t,k,t). (B3) 
e 

It is easy to check that satisfies the same transport equation as where the drift 
term VkW has been replaced by 'S/\^lo — U. 

We now derive the limit of Oe as e ^ 0. Consider the classical asymptotic expansion 

ae~a,a + eai + €^0.2 + . . . . (B 4) 

Upon substitution into (B2) and equating like powers of e, we obtain at order e^^, for 
fixed frequency Wq, 

/ Q(q, k)(ao(x, q, t) - ao(x, k, t))5{uj{ci) - iOo) = 0. (B 5) 

Jq 

It follows from the Krein-Rutman theory (Dautray & Lions (1993) ) that ag is independent 
of q. At order e~^, we obtain 

(Vkco(k) - U) • Vxflo - / Q(q,k)(ai(x,q,t)-ai(x,k,t))5(a;(q)-LJo). (B 6) 

Jq 

The compatibility condition for this equation to admit a solution requires both sides to 
vanish upon integration over d{uj(\i) — ujQ)d\i. Therefore, U satisfies 

V = [ Vku(k)S{to{k)-LUo), where V^, ^ [ S{uj{k) - uo). (B7) 

Once the constraint is satisfied, we deduce from Krein-Rutman theory the existence of a 
vector- valued mean zero corrector x solving 

{Vkc.(k)-U)= / Q(q,k)(x(q)-x(k))<5(c.(q)-c.o) = /:x. (B8) 

Jq 

There is no general analytic expression for which must in practice be solved numer- 
ically. This is typical of problems where the domain of integration in q does not have 
enough symmetries (cf. Allaire & Bal (1999), Bal (1999)). We now have 

ai(x,k,i) =x(k) •Vxao(x,i)- (B9) 

It remains to consider 0(e*^) in the asymptotic expansion. This yields 

ao + (Vkw(k) - tj) • VxSi = £(a2). (B 10) 

The compatibility condition, obtained by integrating both sides over k, yields the wave 
action diffusion equation 

&o-Vx-D-Vx5o = 0, (Bll) 
where the diffusion tensor is given by 

^ = -t!- [ (Vkc.(k) ~ U)x^(k) =-T^ f Cix)x^ik). (B 12) 

The second form shows that D is positive definite since £ is a nonpositive operator. The 



formal asymptotic expansion can be justified rigorously using the techniques in Dautray 
fc Lions (1993)1 . As e ^ 0, we obtain that the error between and oq is at most of 
order e. Therefore, we have that converges to qq satisfying the following drift-diffusion 
equation 

ao + — • Vxao - Vx • D • VxCq = 0, (B 13) 

e 
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with suitable initial conditions. Equation (B 13) is the coordinate-scaled version of (5.16). 
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